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Abstract 

Fungal communities play a key role in ecosystem functioning. However, only 
little is known about their composition in plant roots and the soil of biomass 
plantations. The goal of this study was to analyze fungal biodiversity in their 
belowground habitats and to gain information on the strategies by which ecto- 
mycorrhizal (ECM) fungi form colonies. In a 2-year-old plantation, fungal 
communities in the soil and roots of three different poplar genotypes (Popu- 
lus x canescens, wildtype and two transgenic lines with suppressed cinnamyl 
alcohol dehydrogenase activity) were analyzed by 454 pyrosequencing targeting 
the rDNA internal transcribed spacer 1 (ITS) region. The results were compared 
with the dynamics of the root-associated ECM community studied by morpho- 
typing/Sanger sequencing in two subsequent years. Fungal species and family 
richness in the soil were surprisingly high in this simple plantation ecosystem, 
with 5944 operational taxonomic units (OTUs) and 186 described fungal fami- 
lies. These findings indicate the importance that fungal species are already avail- 
able for colonization of plant roots (2399 OTUs and 115 families). The 
transgenic modification of poplar plants had no influence on fungal root or soil 
communities. Fungal families and OTUs were more evenly distributed in the 
soil than in roots, probably as a result of soil plowing before the establishment 
of the plantation. Saprophytic, pathogenic, and endophytic fungi were the dom- 
inating groups in soil, whereas ECMs were dominant in roots (87%). Arbuscu- 
lar mycorrhizal diversity was higher in soil than in roots. Species richness of 
the root-associated ECM community, which was low compared with ECM 
fungi detected by 454 analyses, increased after 1 year. This increase was mainly 
caused by ECM fungal species already traced in the preceding year in roots. 
This result supports the priority concept that ECMs present on roots have a 
competitive advantage over soil-localized ECM fungi. 



Introduction 

Anthropogenic activities can cause dramatic changes in 
ecosystem structures and their ecological services 
(Dawson 2011). Stability and maintenance of ecosystems 



rely on biodiversity and functional dynamics of organisms 
(Johnson et al. 1996). The impact of organismal groups 
on ecosystem stability depends on several factors such as 
adaptation strategies, interaction with other organisms 
(Johnson et al. 1996), and manner of nutrient acquisition. 
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Fungal Communities of a Poplar Plantation 

Fungi are a group of central importance as they play key 
roles in the carbon and nitrogen cycle improving the 
availability of nutrients for other organisms. They are dis- 
tributed across all climatic zones, and colonize different 
habitats in ecosystems such as soil (Bridge and Spooner 
2001), plant tissues (Arnold et al. 2000), water (Jones 
2011), or rocks (Gadd 2007). 

According to their lifestyle and ecological function, 
fungi can be classified to be saprophytic, pathogenic, 
endophytic, and mycorrhizal. Traditionally, those differ- 
ent groups have been analyzed separately by targeted 
approaches. With the advent of deep sequencing tech- 
niques, it is now possible to record these communities 
comprehensively as a precondition to understanding their 
interactions. For example, the analysis of rhizosphere and 
root endophyte communities in two natural poplar stands 
on contrasting soils revealed differentiation of the com- 
munities between roots and soil as habitats, but surpris- 
ingly, no significant soil-related effects (Gottel et al. 
2011). Furthermore, in contrast to previous morphotyp- 
ing/cloning studies in poplar plantations (Kaldorf et al. 
2004; Stefani et al. 2009), deep sequencing suggested that 
mycorrhiza-forming fungal genera were underrepresented 
in roots (Gottel et al. 2011). It has been speculated that 
genetic differences between poplar species affect mycor- 
rhizal colonization (Tagu et al. 2001; Karlinski et al. 

2010) , and thus influence the composition of fungal com- 
munities in roots (Gottel et al. 2011). Strong variation 
has been found among ectomycorrhizal (ECM) fungi that 
colonize specific coniferous species influenced by plant 
genotypes (Ducic et al. 2009; Karlinski et al. 2010). Many 
ECM fungi show strong host preferences (Lang et al. 

2011) , but the whole root-inhabiting fungal community 
is composed of different ecological groups. It is unknown 
whether fungal root communities as a whole can also be 
affected by the plant genotype. 

Poplars are an important feedstock for biofuel produc- 
tion (Polle and Douglas 2010). Agro-forest areas are cur- 
rently being expanded to meet the demand for 
sustainable biomass production. As soil-borne fungi have 
critical impact on plant health and productivity, the con- 
servation of healthy communities of soil biota and biolog- 
ical soil management are considered pivotal to ensure soil 
fertility and overall productive and sustainable agricul- 
tural systems (Matson et al. 1997). However, knowledge 
of structure, function, and ecology of soil microbial com- 
munities is still very limited, especially for managed agro- 
forest plantations. As there is increasing interest in the 
use of fast growing tree species for production of second 
generation biofuel, attempts are underway to increase pul- 
ping properties of the wood by transgenic modification of 
lignin content and composition (Baucher et al. 1996; 
Pilate et al. 2002). Previous studies showed faster decom- 
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position of leaf litter of poplars with suppressed activity 
of cinnamyl alcohol dehydrogenase (antisense CAD) than 
that of wildtype leaves (Pilate et al. 2002). It is currently 
unknown whether changes in tissue composition of trans- 
genic poplar also influence the assemblage of root-inhab- 
iting fungi or whether transgenic poplars affect the fungal 
community in the soil. 

The main goal of the present study was a compre- 
hensive analysis of fungal biota in soil and roots of 
wildtype and two antisense CAD poplar genotypes to 
test the hypothesis that the soil forms a large species- 
rich reservoir that leads to the differentiation of distinct 
fungal communities in wildtype and transgenic poplars. 
We conducted our study in a recently established exper- 
imental short rotation plantation of hybrid poplar {Pop- 
ulus tremula x P. alba, syn. P. x canescens) wildtype 
and transgenic lines. We applied 454 pyrosequencing 
analyses for in-depth characterization of fungal commu- 
nities using the rDNA ITS1 region as marker gene. The 
role of soil as reservoir for root colonization was inves- 
tigated (i) on the base of taxa composition in fungal 
soil and root communities, (ii) with respect to cluster- 
ing of functional fungal groups in roots of different 
genotypes and adjacent soil, and (iii) with regard to 
temporal dynamics of ECM communities identified by 
morphotyping/sequencing techniques compared to 454 
pyrosequencing. 

Materials and Methods 

Plant materials and study site 

Populus tremula x Populus alba (female clone INRA 
#717-1B4) wildtype and transgenic lines with a modified 
lignin metabolism were multiplied by micropropagation 
(Leple et al. 1992). In June 2008, rooted plantlets were 
planted outdoors in a field trial (47°83'N, 1°91'E) nearby 
the INRA in Orleans, France, on sandy soil with flint. 
Climate is typical of the Loire Valley with oceanic tenden- 
cies, westerly dominant winds, average annual precipita- 
tion of 600 mm, and a mean annual temperature of 10.4° 
C. Natural flora is acidophilic and characteristic of poor 
soils, with oak, birch, chestnut, pine, and heather as 
prominent species belonging to the phytosociologic order 
Quercetalia robori-petraeae. 

The field trial was established in an area of 1365 m 2 
with 120 plants per line (seven transgenic and one wild- 
type) (Fig. 1A). The plants were placed in randomized 
subplots, each consisting of 24 plants (four lines of six 
individual plants) (Supporting Information Fig. SI). The 
plants were drip irrigated during the growing period. In 
March 2010, all trees were coppiced according to typical 
management practices in a short rotation plantation. 
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Fungal Communities of a Poplar Plantation 



Figure 1. Soil cores (B) were taken on a 
2-year-old poplar plantation (A) and cut into 
two longitudinal sections. 454 Pyrosequencing 
was applied on one half to study fungal soil 
and root communities. Out of the second half, 
poplar roots were isolated. Ectomycorrhizal 
fungi colonizing poplar roots were described 
by morphotyping (C) and ITS-sequencing (left 
picture: Hebeloma sacchariolens; right: Laccaria 
tortilis). 




Sampling strategy 

In October 2009, wildtype plants and two transgenic lines 
(ASCAD21 = L21, ASCAD52 = L18) with a decreased 
activity of CAD (Lapierre et al. 1999) were used for sam- 
pling. Three plots per line were chosen (Supporting infor- 
mation Fig. SI). In each plot, nine soil cores (depth: 
0.2 m, diameter: 0.05 m) (Fig. IB) were collected at a 
distance of 0.22 m between two neighboring poplar stems 
(for details, see Supporting Information Fig. SI). In total, 
we collected 81 soil cores (27 per poplar line). In addi- 
tion, leaves were collected. Soil cores and leaves were 
transported on ice and processed in the laboratory within 
72 h after sampling. 

The soil cores (Fig. IB) were cut longitudinally into 
two halves with a sterile scalpel (Supporting information 



Fig. SI). One half was used for analyses of ECM fungal 
community by morphotyping/ITS-sequencing and the 
other half for analyses of the overall fungal soil and root 
community by deep sequencing. 

For ECM analyses three halves were pooled, resulting 
in three samples per plot. The samples were soaked in tap 
water and roots were removed by gentle washing. They 
were stored between wet filter papers at 4°C until further 
processing. 

For analyses of 454 pyrosequencing, each sample was 
processed individually. Roots were cautiously removed 
from the soil, washed in autoclaved water, separated from 
roots of other plant species by shape and color under a 
stereomicroscope (Stemi SV 11, Zeiss, Jena, Germany) 
and frozen at -20°C. The soil was sieved, homogenized, 
subsampled in volumes of 2 mL, and stored at -20°C. 
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Aliquots of the soil samples were used for nutrient ele- 
ment analyses. 

Soil analyses 

Soil pH was determined after extraction in water for 4 h. 
Aliquots of the soil were weighed, dried for 4 days at 60° 
C, weighed again, and used to calculate the dry-to-fresh 
mass ratio. Carbon (C) and nitrogen (N) concentrations 
were determined by dry combustion using a C/N analyzer 
(Carlo Erbas Instruments, Italy). Mineral element concen- 
trations of P, S, K, Ca, Mg, Mn, and Fe were determined 
using an Inductively Coupled Plasma - Atomic Emission 
Spectrometer (Spectro Flame, Spectro Analytic Instru- 
ments, Kleve, Germany) after pressure digestion of sam- 
ples in 65% HN0 3 for 12 h (Heinrichs et al. 1986). To 
determine the nitrate and ammonium concentrations, 
samples of 20 g soil were extracted in 40 ml 1 mmol/L 
CaCl 2 , filtered, freeze-dried, and dissolved in 0.5 mL dou- 
ble deionized water. The aliquots were used for photo- 
spectrometric measurement of nitrate and ammonium 
using commercial kits (Spectroquant, Merck, Darmstadt, 
Germany). 

DNA extraction and quality check 

Eighty-one root samples and 10 leaf samples were freeze- 
dried and ground in a ball mill Type MM2 (Retsch, 
Haan, Germany). Hundred milligrams of root powder 
was suspended in 400 jiL LSS-buffer of the "innuPREP 
Plant DNA kit" (analytikjena, lena, Germany). Genomic 
DNA was extracted according to the manufacturer's 
instructions and eluted in 100 jiL nuclease-free water 
(AppliChem, Darmstadt, Germany). Samples were 
checked for contamination by roots of other plant species 
by amplifying the trnL intron-region of the chloroplast 
DNA with the plant specific primer pair c (CGAA 
ATCGGT AG ACGCT ACG ) and d (GGGGATAGAGG- 
GACTTGAAC) (Taberlet et al. 1991). The polymerase 
chain reaction (PCR) reaction mix was composed of 



2 fiL template DNA (up to 15 ng), 2.5 jiL 10 x buffer 
(Fermentas, St. Leon-Rot, Germany), 2 jiL of MgCl 2 
(25 mmol/L, Fermentas), 1.25 jiL of each primer 
(10 mmol/L) (Eurofms MWG Operon, Ebersberg, Ger- 
many), 0.5 /.(L dNTPs mix (10 mmol/L, Fermentas), 
0.125 /iL Taq polymerase (>10 UlyL, Fermentas), and 
16.625 fiL of nuclease-free water, resulting in a total vol- 
ume of 25 /iL. The PCR was performed in a Mastercycler 
Gradient (Eppendorf, Hamburg, Germany) starting with a 
hot-start at 95°C followed by 95°C for 1 min, 35 cycles of 
30 sec at 94°C (denaturation), 30 sec at 53°C (annealing), 
and 1 min at 72°C (extension), and terminated with 72°C 
for 5 min. PCR products were subjected to electrophore- 
sis in 2% agarose gels, ethidium bromide staining and 
were scanned (Raytest scanner FLA 5100, Straubenhardt, 
Germany). PCR products on the DNA of leaves of the 
same poplar lines as for roots were used as positive con- 
trol. In the few cases where contamination was detected, 
new samples were prepared. 

Eighty-one soil samples were dried in a SpeedVac-Con- 
centrator Savant SPD 11V230 (Thermo, Bonn, Germany) 
and ground in a ball-mill. Genomic DNA was extracted 
using the Soil kit (MoBio, Carlsbad, CA) following the 
manufacturer's instructions. 

Amplicon generation and 454 
pyrosequencing 

All 162 DNA samples were amplified separately. Total 
extracted DNA was employed in the amplification at dif- 
ferent concentrations (undiluted: 1:10, 1:50, 1:100). The 
Amplicon libraries were generated with primers including 
the Roche GS FLX Titanium Amplicon-Adaptor 
Sequences (A-Key, B-Key, Key: TCAG), a 10 bp multiplex 
identifier (MID 1-29, see Table 1, TCB No. 005-2009, 
Roche, Mannheim, Germany) in front of the B-Adaptor 
for multiplexing the PCR products and the template-spe- 
cific primers ITS If (Gardes and Bruns 1993) and ITS2 
(White et al. 1990), resulting in fusion primers A-ITS1F 
(5' CGTATCGCCTCCCTCGCGCCATCAG-CTTGGTCAT 



Table 1. Summary of 454 pyrosequencing data. 





Soil beneath 




Roots from 






Transgenic poplar 


Wildtype poplar 


Transgenic poplar 


Wildtype poplars 


Sequence reads 

Sequence reads after quality control 
Sequences/sample 

Number of OTUs (nonsingletons/sample) 
Number of singletons/sample 


297,836 
251,883 
11,631-15,965 
392-800 
326-675 


153,626 
129,962 
9524-17,994 
395-736 
307-703 


203,238 
166,556 
6568-10,835 
75-225 
48-143 


157,200 
137,652 
4706-9620 
1 1 8-249 
72-112 



Samples are defined by sample type (soil or root samples) and poplar genotype (transgenic or wildtype). Twenty-seven samples were taken per 
sample type: Eighteen samples of transgenic and nine of wildtype plants, respectively. OTUs, operational taxonomic units 
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TTAGAGGAAGTAA- 3') and B-MID-ITS2 (5' CTATGC 
GCCTTGCCAGCCCGCTCAG-MID-GCTGCGTTCTTCA 
TCGATGC). PCR reactions were performed as described 
above, but 0.7 /iL of 16 mg/mL bovine serum albumin 
(Merck, Darmstadt, Germany) was added to a total PCR 
mix volume of 25 fiL. After amplification, the PCR prod- 
ucts were purified using the "innuPREP PCRpure Kit" 
(analytikjena). Then, the PCR products from three cores 
of the same tree were pooled, resulting in 27 Amplicon 
libraries for root and soil, respectively, with independent 
replicates for each line. Amplicon concentration was 
determined using the Qubit M dsDNA HS Assay Kit in a 
Qubit fiuorometer (Invitrogen GmbH, Karlsruhe, Ger- 
many). The 27 amplicon libraries of root- and soil-sam- 
ples, respectively, were pooled in equal amounts for 454 
pyrosequencing. Amplicon libraries were sequenced with 
the 454 Genome Sequencer FLX (Roche) using the ampli- 
con sequencing protocol and Titanium chemistry 
(Roche). Sequencing was performed by the Gottingen 
Genomics laboratory (http://www.g2l.bio. uni-goettingen. 
de/). Three medium lanes of a Titanium pico titer plate 
were used for sequencing of the complete amplicon 
libraries. The entire quality data set as unprocessed data 
files were deposited into the sequence read archive (SRA). 
The study accession number is ERP001442 and can be 
accessed by following link: http://www.ebi.ac.uk/ena/data/ 
view/ERP001442. 

Bioinformatics and OTU clustering 

After the removal of barcodes and tags, 454 pyro sequenc- 
ing reads were processed with a perl script discarding all 
reads shorter than 150 bp and reads containing more 
than four ambiguity symbols. On average, 74% of all 
reads passed these criteria. The individual sample FASTA 
files were subjected to cluster analysis for a tentative OTU 
count using the clustering function of USEARCH v5.3.23 
(Edgar 2010) with the following criteria: > 97% similarity 
over > 90% sequence length. Cluster analyses were car- 
ried out on individual and "combined" samples, the latter 
ones including sequence read information of all soil or all 
root samples. All singletons were removed prior to further 
analyses. To identify OTUs at taxonomic level, a ran- 
domly selected sequence of an OTU cluster was compared 
with the nonredundant GenBank database (Benson et al. 
2008) and the custom-curated database RSyst (http://my- 
cor.nancy.inra.fr/RSyst/) using BLASTn (Altschul et al. 
1997). A postprocessing perl script stored the 10 best 
BLASTn hits per cluster with an expectation value of 
<10e~ 3 in a BLASTn-file. OTUs with a taxonomic assign- 
ment at the species level were classified with respect to 
their ecological lifestyle by literature research (Table S3, 
Supplemental information). Ecological groups were cate- 
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gorized as follows: AM, arbuscular mycorrhizal; ECM, 
ectomycorrhizal; lichenized; saprotrophic; endophytic or 
parasitic. 

Morphotyping on root tips 

Grass roots were identified by differences in morphology 
and removed. Three-hundred living root tips were 
inspected per poplar tree. ECM fungi were morphotyped 
(Fig. 1C) using a simplified method after Agerer (1987- 
2006) recording shape, color, texture of the mantle, and 
presence or absence of hyphae or rhizomorphes under a 
stereomicroscope (M205 FA, Leica, Wetzlar, Germany). 
ECM colonization (%) was calculated as: number of ECM 
root tips x 100/total number of root tips. Three to four 
ECM root tips of each morphotype were collected and 
stored at -20°C. 

Cloning and sequencing of ECM species 

Genomic DNA of the frozen ECM root tips was 
extracted using the "innuPREP Plant DNA kit" (analy- 
tikjena). The rDNA ITS-region was amplified by PCR 
with the primer pair ITS5/ITS4 (White et al. 1990) as 
described above with the following modifications: 34 
cycles and an annealing temperature of 55°C. Direct 
Sanger-sequencing or cloning/sequencing was carried out 
according to Lang et al. (2011). Sequences were blasted 
using the following databases: NCBI (nBLAST) (http:// 
www.ncbi.nih.gov/), Fungal RSyst (http://mycor.nancy. 
inra.fr/RSyst/), and UNITE (http://unite.ut.ee/). 
Sequences are available at NCBI (accession JQ409279- 
JQ409296). 

Data analyses 

To test for possible variability of fungal communities of 
different samples, 454 pyrosequencing data were blasted 
against the RSyst database. A perl script stored the top 
BLASTn hit (E-value < 10e~ 3 ) and the number of reads 
per species of each sample in a csv-file. Statistical analy- 
ses were performed on the basis of the number of reads 
per species and the relative abundance of reads. Samples 
were compared by a pairwise test based on the relative 
frequencies. The Wilcoxon rank-sum test was used 
to identify significant differences according to a 
P-value < 0.05 after Bonferroni-correction. All our 
statistical analyses were carried out by using the software 
R-2.9.2 (R Development Core Team 2009). Additionally, 
nonmetric multidimensional scaling (NMDS) with the 
function metaMDS of the "vegan" package was applied. 
Before running the NMDS, data were square root trans- 
formed. 
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Statistical analyses of the fungal communities forming 
visible ECMs with roots were based on the relative abun- 
dance of the morphotypes. The Kruskal-Wallis rank sum 
test (package "stats") was used to identify differences 
between poplar genotpyes (P < 0.05). 

The Wilcoxon rank sum test with an additional Bonfer- 
roni correction was carried out to examine differences in 
biodiversity indices and the relative abundances of fungal 
families within different sample types, respectively. 

The defined OTUs were used to calculate taxon accu- 
mulation curves with the freeware software Analytic Rare- 
faction version 1.3 (http://www.uga.edu/strata/software/ 
Software.html). Biodiversity indices and species richness 
estimators were calculated using the software Estimates 
version 8.0.0 (Collwell 2006). Evenness was additionally 
determined by the formula (Shannon/LN [number of 
detected OTUs]). 

Presence/absence data of fungal families in individual 
soil and root samples were subjected to hierarchical clus- 
ter analysis using EPCLUST (http://www.bioinf.ebc.ee/EP/ 
EP/EPCLUST/index.cgi). Correlation-based distance mea- 
sure was chosen as similarity metric and average distance 
as clustering method. 

Differences in soil parameters were tested with one-way 
analysis of variance ( ANO VA) . 

Results 

Fungal species richness and diversity in soil 
and roots 

In total, 811,900 sequence reads were generated by 454 
pyrosequencing. Sequence reads that did not match our 
quality criteria were removed (see Material and Methods) 
resulting in 686,053 sequence reads for further analyses. 
In all, 4706-17,994 sequences were obtained per sample 
(Table 1). These sequences were clustered according to 
similarity and yielded 75-800 nonsingleton OTUs per 
sample (Table 1). Forty-eight to 703 singletons per sam- 
ple were obtained (Table 1). 

Rarefaction curves based on 97% sequence identity lev- 
eled off between 398 and 817 OTUs for soil samples and 
between 91 and 249 OTUs for root samples (Supporting 
Information, Fig. S2A and B). Rarefaction analyses for 
complete fungal richness of the study site in soil and 
roots showed saturation at 5944 and 2399 OTUs, respec- 
tively (Supporting Information, Fig. S2D). Root samples 
exhibited higher variability in the shape of their species 
accumulation curves indicating strong scattering of spe- 
cies richness between different samples. Rarefaction analy- 
ses of ECM root communities revealed complete coverage 
(Supporting Information, Fig. S2C). Estimated species 
richness (H max ) showed a clear decrease in the order of 




H(max) H' E 



Figure 2. Species richness (H max ), Shannon index (H'), and Evenness 
(E) of fungi in soil, roots, and of root-associated ECM communities. 
Diversity indices are means (n = 27 ± SE). Significant differences 
between bar heights (P < 0.0S) are indicated by different letters 
above bars. H max = In (species number). 

the habitats soil > roots > root-associated ECM commu- 
nities (Fig. 2). The decrease in species richness from soil 
to roots and ECM communities was also reflected by the 
Shannon indices (Fig. 2). Evenness was highest for ECM 
communities and lowest for fungi in roots (Fig. 2). 

Fungal community structure in different 
habitats and poplar genotypes 

To find out whether the poplar genotype affected fungal 
abundance or community structures in roots or soil, Wil- 
coxon rank-sum tests with Bonferroni correction were 
conducted. However, no significant difference was 
detected among the three investigated genotypes (trans- 
genic lines ASCAD52 [=L18], ASCAD21 [=L21] and wild- 
type; P < 0.05) with respect to the presence of fungal 
species or their abundance. These findings held true for 
soil and root samples as well as root-associated ECM 
fungi. Soil nutrients (per gram dry soil) did not differ 
between samples of different poplar genotypes (N0 3 ~, 
17.2 ± 2.04 /imol; NH 4 + , 15.3 ± 1.7 junol; total N, 0.88 ± 
0.08 mg; P, 0.22 ± 0.01 mg; S, 0.09 ± 0.01; Ca, 0.93 ± 
0.04 mg; Mg, 0.51 ± 0.02 mg; Mn, 0.17 ± 0.01 mg; Fe, 
4.16 ± 0.20 mg; C, 15.5 ± 1.4 mg; pH 5.85 ± 0.03), with 
the exception of K (mean: 1.19 ± 0.04 mg/g dry soil), 
which was slightly higher (13% above the mean) in soil 
collected beneath poplar line 18 than in that beneath the 
wildtype (P < 0.02). 

An NMDS plot calculated for OTUs revealed strong 
clustering of fungal communities for soil and roots, 
respectively (stress = 13.63, nonmetric fit R 2 = 0.98) 
(Fig. 3). Permutation test confirmed significant classifica- 
tion with P < 0.001 (R 2 = 0.6332). No separation of 
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Figure 3. NMDS plot of the fungal community structure using the 
Bray-Curtis dissimilarity measure. Each point represents the fungal 
community of a given sample. Permutation tests revealed a highly 
significant classification (P = 0.001). Samples were classified according 
to the plant genotype (wildtype; L18 and L21, transgenic), sampling 
point, and sample type (black circles, soil; open circles, root). Stress 
value =13.63, R 2 = 0.98. Ellipses separate samples into two 
categories: left ellipse, soil samples; right one, root samples. 
Confidence area of ellipses = 0.95. 

samples related to plant genotype or the position in the 
field was detected (Fig. 3). The significant differences 
between the fungal communities of soil and roots origi- 
nated, therefore, from the lower species richness of roots 
compared with soil. Although soil contained higher spe- 
cies richness than roots, the scattering of data was lower 
(see ellipses in Fig. 3), indicating higher homogeneity of 
species distribution in soil than in roots. 

Fungal family abundance and distribution 
across soil and root samples 

OTUs were clustered according to their taxonomic affilia- 
tion into overall 196 fungal families. Soil (186) and root 
(115) samples differed in their fungal family composition. 
Eighty-one fungal families were solely found in soil sam- 
ples and 10 only in root samples (Supporting Informa- 
tion, Table S2). Among the common families, 59 were 
significantly more abundant in soil samples than in roots 
(Fig. 4). Two families, Filobasidiaceae and Mortierellaceae, 
were dominant in soil, each comprising about 15% of all 
OTUs. The relative abundance of nine further families in 
soil ranged between 1% and 5%, whereas all other fungal 
families that differed significantly from roots were present 
only with low abundance (<1%) (Fig. 4). 

In root samples, six fungal families were significantly 
enriched in comparison with soil (Fig. S3). Pyronemata- 
ceae dominated the community (13.5%) in roots, while 
the relative abundances of Paxillaceae, Paraglomeraceae, 
Rhytismataceae, and Sporormiaceae ranged between 1.3% 




0.0 0.5 1.0 1.5 2.0 2.5 10 20 30 40 

Relative abundance in soil (%) 



Figure 4. Distribution of fungal families in soil and root samples 
according to their relative abundance. Significant differences 
(P < 0.05) between soil and roots are indicated by black circles; open 
circles indicate fungal families with similar abundance in soil and 
roots. Families with abundances above 1 % were labeled with a star. 
Red line indicates equal abundances in both roots and soil; 100% is 
the total abundance of all fungal families. 



and 2%. Russulaceae were represented by 0.04% of the 
OTUs (Supporting Information, Table SI). 

Hierarchical cluster analyses demonstrated the distribu- 
tion pattern of fungal families in individual samples 
(Fig. 5). In soil samples, about one quarter of all fungal 
families were present in >90% of the samples. Forty-six 
percent and 59% of fungal families were detected in at 
least >50% and >25% of all soil samples, respectively 
(Fig. 5A). In contrast, the clustering of fungal families in 
root samples differed (Fig. 5B). Only 8% of all fungal 
families in root samples were present in >90% of all sam- 
ples. Twenty percent and 38% of fungal families were 
present in >50% and >25% of the samples, respectively 
(Fig. 5B). 

Ecological groups in soil and root samples 

To examine the distribution of ecological groups in soil 
and root samples 1272 and 463 OTUs, respectively, that 
could be assigned to species levels were selected and their 
abundances were set to 100%. One hundred and fifty-six 
and 27 of these species constituted 90% of the relative 
abundance in soil and roots, respectively, and were classi- 
fied after literature research as ECM, AM, saprophytic, 
endophytic, pathogenic, or lichenized fungi (Supporting 
Information, Table S3). In soil samples, saprophytic fungi 
(47%) formed the largest group, followed by 23% ECM, 
19% pathogenic, and 8% endophytic fungi. Lichenized 
and AM fungi were present only in low abundances of 
1.8% and 0.4%, respectively (Fig. 6). 
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(A) (B) 




Figure 5. Heat map showing clustering of fungal families in (A) soil and (B) root samples. The color code of the heat map indicates presence 
(red) or absence (black) of fungal families (in rows) in the individual samples (in columns). 



In root samples, ECM fungi were the dominant 
group encompassing 87% of the total abundance. Endo- 
phytic, parasitic, and saprophytic fungi showed lower 



abundances of 5%, 4%, and 4%, respectively. On the 
species level, no AM or lichenized fungi were detected 
(Fig. 6). 
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Figure 6. Distribution of fungal species with different ecological lifestyle in soil and root samples. OTUs with >97% sequence identity to known 
species were classified based on their taxonomic affiliation to six different ecological lifestyles. Only OTUs accounting for 90% of overall relative 
abundance were included in the analysis. All ecological lifestyle groups were significantly different (P < 0.05) from the dashed line, which 
indicates the same abundance in roots and soil. AM, arbuscular mycorrhizal; ECM, ectomycorrhizal fungal species; OTU, operational taxonomic 
unit. 



Dynamic of the ECM community on poplar 
roots 

With increasing age, poplar roots showed a typical 
increase in ECM fungal richness (Smith and Read 2008). 
In October 2009, seven, and in October 2010, nineteen 
ECM fungal species were detected on roots, of which six 
(2009) and 16 (2010) were identified by rDNA ITS 
sequencing (Supplemental information Table S4). With 
the exception of Hebeloma sp. and an uncultured Peziz- 
ales (DQ469743.1), the ECMs identified in 2009 were also 
present in 2010. 

To understand dynamic processes within the ECM 
community and root colonization, morphotyping/ITS- 
sequencing and 454 pyrosequencing approaches were 
compared. All ECM species detected in 2009 were also 
detected by 454 pyrosequencing in both soil and root 
samples (Table 2). Furthermore, 13 of the 16 fungal spe- 
cies that colonized the roots in 2010 were already detected 
on poplar roots by 454 pyrosequencing in 2009. Two of 
the three missing species, Scleroderma bovista, and Tuber 
sp., were detected solely in soil samples. Only one species, 
an uncultured Ascomycota (EU557319.1) that formed an 
ECM in 2010, was detected neither in soil nor root sam- 
ples in 2009 by 454 pyrosequencing approach. 



Discussion 

Massive 454 pyrosequencing reveals 
surprisingly high fungal species richness in 
a young short rotation plantation 

Rarefaction analyses indicated that we detected the major- 
ity of nonsingleton OTUs present in soil (average 556 per 
sample, 5944 OTUs for the complete survey) and roots 
(145 per sample, 2399 OTUs) of the complete experimen- 
tal site of a 2-year-old poplar stand (Supporting Informa- 
tion, Fig. S2A and B). These numbers are relatively high 
compared with other studies reporting deep sequencing 
of fungal communities in soil of mature forest stands 
(Buee et al. 2009; Gottel et al. 2011) and roots from 
mature oak trees (Jumpponen et al. 2010). One reason 
may be a higher sampling density in our study compared 
with the previous ones. Nevertheless, it is remarkable that 
even in simple and young agro-ecosystems established on 
a tilled soil (Fig. 1A), very high sequencing depth is 
needed for comprehensive characterization of fungal com- 
munity composition. 

The fungal family richness (186 in soil and 115 in pop- 
lar roots) also exceeded values that have been previously 
reported for fungal soil communities (O'Brien et al. 2005; 
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Table 2. Fungal species detected by two approaches: morphotyping/ 
ITS-sequencing and 454 pyrosequencing. 





ITS-sequencing 


454-Pyrosequencing 




Ectomycorrhizal 








poplar 


Poplar 






root tips 


Soil 




roots 


Fungal species 


2009 2010 


2009 


2009 


Peziza ostracoderma 


X X 






Paxillus involutus 


X X 


X 


X 


Laccaria tortilis 


X X 






Hebeloma 


X X 


X 


X 


sacchariolens 








Tomentella ellisii 


X 






Scleroderma bovista 


X 




X 


Cenococcum geophilum 


X 






Xerocomus ripariellus 


X 


X 


X 


Hebeloma sp. 


X 


X 


X 


Geopora sp. 








Tuber sp. 


X 




X 


Uncultured Ascomycota 


X 






JQ409293 








Uncultured Ascomycota 


X 


X 


X 


JQ409292 








Uncultured 


X 


X 


X 


ectomycorrhizal 








funni ID40Q7Q4 

uiiyi J^^+UIJZ. -?^+ 








Uncultured fungus 


X 


X 


X 


JQ409288 








Uncultured fungus 


X 


X 


X 


JQ409287 








Uncultured Peziza 


X 


X 


X 


JQ40929B 








Uncultured Pezizales 


X 


X 


X 


JQ409284 









In October 2009 and 2010, ectomycorrhizal poplar root tips were 
sampled, classified by morphotyping and analyzed by ITS-sequencing 
(in total 27 samples). Additionally in 2009, poplar roots and soil sam- 
ples were taken and subjected to 454 pyrosequencing analysis. 



Buee et al. 2009), fungal phyllosphere (Jumpponen and 
Jones 2009), and root communities of oak (Jumpponen 
et al. 2010) (Supporting Information, Table SI). As no 
adjacent forest or mature site existed that could cause 
"vicinal invasion" (Kaldorf et al. 2004), our study shows 
that already very young stands own a rich and diverse 
reservoir of fungal propagules. 

Roots and soil constitute distinct ecological 
fungal biomes 

We observed a clear separation of soil and root fungal 
communities (Fig. 3). A clear separation of microbiomes 
has also been reported for the rhizosphere and endo- 
sphere of mature poplar sites (Gottel et al. 2011). Our 
study shows that the differentiation of these habitats 



occurs already in an early phase of stand development 
and is mainly the result of fungal families enriched in soil 
(about 1/3 of all soil families) compared with roots. This 
observation points to high selectivity of interactions of 
roots with soil fungal genera (Fig. 3; Supporting Informa- 
tion, Table SI). The majority of significant fungal soil 
families shared saprophytic or parasitic lifestyles (Fig. 6) 
including the two most abundant fungal soil families, the 
Filobasidiaceae and Mortierellaceae (Hibbett et al. 2007). 
Members of these families are widespread, occurred also 
with high abundance in soils of six different tree mono- 
plantations and have therefore been classified as generalis- 
tic families (Buee et al. 2009). 

Analysis of the lifestyle of the most abundant fungal 
species revealed significant enrichment of pathogenic, 
endophytic, lichenized, and AM fungi in soil compared 
with roots (Fig. 5). Some earlier studies demonstrated 
that pathogenic fungi are forming a large group within 
fungal communities in plant tissues (Monk and Samuels 
1990; Bills and Polishook 1994) and that (bacterial) 
antagonists affect overall abundance of pathogenic fungi 
(Berg et al. 2002). However, the analysis of fungal com- 
munities in plant tissue samples has been challenging in 
the past due to inadequate isolation techniques (Bayman 
2007). For example, in a deep sequencing study, Jumppo- 
nen et al. (2010) reported 12.3% of all detected fungi in 
mycorrhizal oak roots to be pathogenic. Our study shows 
that fungi with this lifestyle were about five-times more 
abundant in soil than in roots (Fig. 5). 

Interestingly, the abundance of endophytic fungi was 
also higher in soil than in roots (Fig. 6). The mechanism 
of endophytic transmission is very variable and depends 
on the endophytic class (Rodriguez et al. 2009) ranging 
from spores dispersed by wind or rain to released hyphal 
fragments or infected (dead) plant tissue passively distrib- 
uted by herbivores (Monk and Samuels 1990) or physical 
disturbance. These pathways and the influence of abiotic 
factors such as land-use lead to sometimes unexpected 
abundance and diversity of endophytic fungi (Rodriguez 
et al. 2009) found in different biomes such as agro-sys- 
tems and terrestrial ecosystems (Arnold and Lutzoni 
2007). Additionally, the identification of fungi as endo- 
phytes is problematic as the classification is often based 
on the momentary status of detection without regarding 
the future status of interaction (Schulz and Boyle 2005). 
Thus, fungi termed endophytic might be saprophytic or 
pathogenic in a certain part of their lifecycle. 

Some distinct classes of mutualistic fungi including two 
families of mycorrhizal fungi (Archaeosporaceae [AM] and 
Bankeraceae [ECM]; Fig. 6, Supporting Information, 
Table SI) were significantly enriched in soil. The over- 
abundance of AM fungi in soil is surprising as poplar trees 
are able to associate with both AM and ECM fungi at the 
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same time (Molina et al. 1992). However, here ECM fungi 
formed the largest ecological group in roots with almost 
90% abundance (Fig. 6) more than previously reported by 
Jumpponen et al. (2010) for ECM-colonized oak roots 
(72%). The strong colonization with ECM was probably 
caused by preceding long-term cultivation of poplars on 
the experimental sites and this may have suppressed AM 
proliferation (Dhillion 1994; Chen et al.2000). 

The ECM accumulation in roots was mainly due to 
OTUs assigned to four families: Inocybaceae, Pezizaceae, 
Paxillaceae, and Pyronemataceae (Fig. 4). Whereas the for- 
mer two were evenly distributed between soil and roots, 
the latter two were predominantly present in roots. 
Assignment of pezizalean Pyronemataceae taxa to specific 
ecological lifestyles remained problematic longtime, as they 
comprise a heterologous family. In fact, they are nowadays 
considered paraphyletic (Perry et al. 2007). In our study, 
Pyronemataceae showed significant presence in root sam- 
ples (Fig. 4) and were one of the families with the highest 
genera richness (Supporting Information, Fig. S3). The 
different distribution of genera in soil and roots supports 
previously assigned ecological lifestyles of some taxa of the 
Pyronemataceae: in roots solely, genera described as 
mycorrhiza forming fungi were detected, while in soil, 
additional taxa with other ecological lifestyles were found. 

The distribution of fungal families in individual sam- 
ples was more homogeneous in soil than in roots (Fig. 5). 
This was also supported by the narrow clustering of 
OTUs in the NMDS analysis (Fig. 3) and the larger calcu- 
lated Evenness in soil than in roots (Fig. 2). Mycorrhizal 
fungi are known to cluster along the root system of their 
host plants forming a patchy distribution (Smith and 
Read 2008). This may also be expected for fungal soil 
communities on early successional sites, as soil factors 
can differ widely at one site (Reverchon et al. 2010). In 
our study, the lack of significant differences in soil factors 
and soil plowing before the establishment of the planta- 
tion may have resulted in the relative homogeneous dis- 
tribution of soil inhabiting fungi. The observation that a 
small number of ECM forming genera were dominant in 
roots and that roots contained a high number of rare 
OTUs at the same time suggests that roots were underly- 
ing high colonization pressure, but fungal proliferation 
was effectively suppressed with the exception of ECM. 
However, further studies are needed to shed light on the 
mechanisms influencing the composition of ecological 
groups in fungal communities in different habitats. 

Deep sequencing reveals host effects on the 
priority of ECM root colonization 

The application of a double approach, morphotyping/San- 
ger-sequencing and 454-pyrosequencing, allowed us to 
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draw a picture of dynamic processes and cross-links of 
fungal soil and root communities in relation to ECM col- 
onization. The ECM community on poplar roots showed 
the well-known increases in colonization rate and diver- 
sity with increasing tree age (Dhillion 1994; Chen et al. 
2000; Egerton-Warburton and Allen 2001). The fungal soil 
community (2009) harbored already all but one of the 
fungal species that formed ECM with poplar roots in 
the following year (2010, Table 2). This finding indicates 
the strength of fungal soil communities as a source for 
plant root colonization, and suggests low invasion by soil 
fungi from outside the agro-system within one annual 
cycle. Furthermore, most fungal species with ECM devel- 
opment in 2010 were already traced on poplar roots in 
2009 by 454 pyrosequencing (Table 2). The experimental 
site was underlying early successional dynamics with fac- 
tors that influence fungal root colonization such as the 
pattern of C allocation (Druebert et al. 2009), fungal 
competition (Kennedy et al. 2009), or availability of 
nutrients (Peter et al. 2001). While pronounced changes 
in soil nutrient availability appear unlikely, the growth of 
the poplars from about 0.2 to 1.9 m in the first year after 
planting (L. Danielsen, unpubl. results) indicates a strong 
increment in carbon productivity, which is one of the 
main drivers of ECM diversity (Druebert et al. 2009; Pena 
et al. 2010). The priority concept for ECM colonization, 
which has experimental support under controlled condi- 
tions (e.g., Kennedy et al. 2009), holds that the first 
mycorrhizal species to colonize a host's roots subsequently 
is the stronger competitor, when other fungal species are 
added. Our results suggest that this concept needs to be 
expanded to account for the dynamics of the colonized 
habitat. Most changes in ECM root communities were 
caused by fungal species already present on roots, that is, 
prior to other ECM present in soil that became more 
competitive forming functional ECM in the second year. 
As there were no changes in climatic or edaphic factors, 
which could have resulted in changes in the ECM assem- 
blages, plant-related factors such as changes in carbon 
availability must have been responsible for the shift in the 
dominance of fungal species in the ECM communities. 

Transgenic poplars with suppressed CAD 
activity do not affect soil, root, or ECM 
communities 

One important goal of this study was the assessment of 
the impact of transgenic versus wildtype poplar plants on 
fungal soil, root, and ECM communities; but no signifi- 
cant differences were observed (Fig. 3). Previous studies 
have already indicated no influence of transgenic poplar 
genotypes (rolC, a transformation causing stunting; npll:: 
GUS, a selection marker coupled with a reporter gene) on 



© 2012 The Authors. Published by Blackwell Publishing Ltd. 



1945 



Fungal Communities of a Poplar Plantation 

ECM community structures (Kaldorf et al. 2004; Stefani 
et al. 2009). Here, we show that this also holds true for 
transgenic poplars (antisense CAD) with improved pul- 
ping properties that were modified in their phenylpropa- 
noid metabolism (Pilate et al. 2002). This is an important 
result because other studies revealed significant correla- 
tions between genetically altered phenolic concentrations 
and associated above-ground organismic interactions 
(Kleemann et al. 2011). Earlier studies on genetically 
modified poplars were limited because only ECM or culti- 
vable soil fungi could be analyzed. Our data add impor- 
tant information with regard to the bio-safety discussion 
because we show that in situ fungal soil and root com- 
munities were unaffected by host modification of an 
important commercial trait. These results are especially 
interesting where fungi are concerned that depend on 
host plant features, such as endophytic or parasitic fungi. 
In contrast to our working hypothesis, we did not detect 
any significant differences between the fungal communi- 
ties of wildtype and antisense CAD poplars. Nevertheless, 
it is clear that genotype x biotic environment interac- 
tions cannot be excluded in general because intraspecific 
variations of ECM colonization have been demonstrated 
in crossing pedigrees (Labbe et al. 2011). Therefore, biotic 
interactions will have to be tested for each transgenic line 
that is planted in the field. 

Conclusions 

The results of our analyses indicate that fungal soil and 
root community interact by dynamic processes and that 
soil is playing an important role as a fungal reservoir. 
Poplar roots were dominated by ECM fungi. The downre- 
gulation of an enzyme of lignin biosynthesis (antisense- 
CAD) did not affect ECM, root, or soil fungal 
assemblages. To our knowledge, we described for the first 
time the proportional composition of fungal ecological 
groups of two interacting fungal communities. Informa- 
tion on ecological groups and composition of fungal 
communities is urgently needed to understand the vari- 
able nature of fungal communities and underlying mecha- 
nisms of interaction. Additionally, the combination of 
two different detection techniques allowed us to draw a 
comprehensive picture of fungal soil and root communi- 
ties of the experimental site. 

Acknowledgments 

We are grateful to M. Reichel and T. Klein (Laboratory 
for Radio-Isotopes) for help with sampling, root harvest 
and DNA extraction, to D. Janz for help with OTU clus- 
tering and to L. Reich for English proofreading. We grate- 
fully acknowledge financial support by the European 



L. Danielsen et al. 

Commission within the Seventh Framework Program for 
the Research, Project Energypoplar (FP7-211917). 

Data Accessibility 

454 Pyrosequencing raw data set: SRA study accession 
number: ERP001442 (http://www.ebi.ac.uk/ena/data/view/ 
ERP001442). 

Conflict of interest 

None declared. 
References 

Altschul, S. F., T. L. Madden, A. A. Schaffer, J. H. Zhang, Z. 

Zhang, W. Miller, et al. 1997. Gapped BLAST and PSI- 

BLAST: a new generation of protein database search 

programs. Nucleic Acids Res. 25:3389-3402. 
Arnold, A. E., and F. Lutzoni. 2007. Diversity and host range 

of foliar fungal endophytes: are tropical leaves biodiversity 

hotspots? Ecology 88:541-549. 
Arnold, A. E., Z. Maynard, G. S. Gilbert, P. D. Coley, and T. 

A. Kursar. 2000. Are tropical fungal endophytes 

hyperdiverse? Ecol. Lett. 3:267-274. 
Baucher, M., B. Chabbert, G. Pilate, J. van Doorsselaere, M. T. 

Tollier, M. Petit-Conil, et al. 1996. Red xylem and higher 

lignin extractability by down-regulating a cinnamyl 

alcohol dehydrogenase in poplar. Plant Physiol. 

112:1479-1490. 

Bayman, P. 2007. Fungal endophytes. Springer-Verlag, Berlin 
Heidelberg. 

Benson, D. A., I. Karsch-Mizrachi, D. J. Lipman, J. Ostell, and 
D. L. Wheeler. 2008. GenBank. Nucleic Acids Res. 
36:D25-D30. 

Berg, G., N. Roskot, A. Steidle, L. Eberl, A. Zock, and K. 

Smalla. 2002. Plant-dependent genotypic and phenotypic 

diversity of antagonistic rhizobacteria isolated from different 

Verticillium host plants. Appl. Environ. Microbiol. 

68:3328-3338. 
Bills, G. F., and J. D. Polishook. 1994. Abundance and 

diversity of microfungi in leaf-litter of a lowland rain-forest 

in Costa-Rica. Mycologia 86:187-198. 
Bridge, P., and B. Spooner. 2001. Soil fungi: diversity and 

detection. Plant Soil 232:147-154. 
Buee, M., M. Reich, C. Murat, E. Morin, R. H. Nilsson, S. 

Uroz, et al. 2009. 454 Pyrosequencing analyses of forest soils 

reveal an unexpectedly high fungal diversity. New Phytol. 

184:449-456. 

Chen, Y. L., M. C. Brundrett, and B. Dell. 2000. Effects of 
ectomycorrhizas and vesicular-arbuscular mycorrhizas, alone 
or in competition, on root colonization and growth of 
Eucalyptus globulus and E. urophylla. New Phytol. 146: 
545-555. 



1946 



© 2012 The Authors. Published by Blackwell Publishing Ltd. 



L. Danielsen et al. 



Fungal Communities of a Poplar Plantation 



Collwell, R. K. 2006. EstimateS: statistical estimation of species 
richness and shared species from samples. Version 8.0. User' 
s guide and application. Available at http://viceroy.eeb. 
uconn.edu/estimates/index.html (accessed June, 21, 2012). 

Dawson, T. P. 2011. Beyond predictions: biodiversity 

conservation in a changing climate (vol 332, pg 53, 2011). 
Science 332:664-664. 

Dhillion, S. S. 1994. Ectomycorrhizae, arbuscular mycorrhizae, 
and Rhizoctonia sp. of alpine and boreal Salix spp. in 
Norway. Arct. Alp. Res. 26:304-307. 

Druebert, C, C. Lang, K. Valtanen, and A. Polle. 2009. Beech 
carbon productivity as driver of ectomycorrhizal abundance 
and diversity. Plant Cell Environ. 32:992-1003. 

Ducic, T., D. Berthold, R. Langenfeld-Heyser, F. Beese, and A. 
Polle. 2009. Mycorrhizal communities in relation to biomass 
production and nutrient use efficiency in two varieties of 
Douglas fir (Pseudotsuga menziesii var. menziesii and var. 
glauca) in different forest soils. Soil Biol. Biochem. 41: 
742-753. 

Edgar, R. C. 2010. Search and clustering orders of magnitude 
faster than BLAST. Bioinformatics 26:2460-2461. 

Egerton-Warburton, L., and M. F. Allen. 2001. Endo- and 
ecotmycorrhizas in Quercus agrifolia Nee. (Fagaceae): 
patterns of root colonization and effects on seedling growth. 
Mycorrhiza 11:283-290. 

Gadd, G. M. 2007. Geomycology: biogeochemical 

transformations of rocks, minerals, metals and radionuclides 
by fungi, bioweathering and bioremediation. Mycol. Res. 
111:3-49. 

Gardes, M., and T. D. Bruns. 1993. ITS primers with enhanced 
specificity for Basidiomycetes - application to the 
identification of mycorrhizae and rusts. Mol. Ecol. 2: 
113-118. 

Gottel, N. R., H. F. Castro, M. Kerley, Z. M. Yang, D. A. 
Pelletier, M. Podar, et al. 2011. Distinct microbial 
communities within the endosphere and rhizosphere of 
populus deltoides roots across contrasting soil types. Appl. 
Environ. Microbiol. 77:5934-5944. 

Heinrichs, R., H. J. Brumsack, N. Loftfield, and N. Konig. 
1986. Improved pressure digestion system for biological and 
inorganic materials. Z. Pflanzenrnahrung. Bodenk. 149: 
350-353. 

Hibbett, D. S., M. Binder, J. F. Bischoff, M. Blackwell, P. F. 
Cannon, O. E. Eriksson, et al. 2007. A higher-level 
phylogenetic classification of the Fungi. Mycol. Res. 
111:509-547. 

Johnson, K. H., K. A. Vogt, H. J. Clark, O. J. Schmitz, and D. 

J. Vogt. 1996. Biodiversity and the productivity and stability 

of ecosystems. Trends Ecol. Evol. 11:372-377. 
Jones, E. B. G. 2011. Fifty years of marine mycology. Fungal 

Diversity 50:73-112. 
Jumpponen, A., and K. L. Jones. 2009. Massively parallel 454 

sequencing indicates hyperdiverse fungal communities in 



temperate Quercus macrocarpa phyllosphere. New Phytol. 
184:438-448. 

Jumpponen, A., K. L. Jones, J. D. Mattox, and C. Yaege. 2010. 
Massively parallel 454-sequencing of fungal communities in 
Quercus spp. ectomycorrhizas indicates seasonal dynamics 
in urban and rural sites. Mol. Ecol. 19:41-53. 

Kaldorf, M., C. Renker, M. Fladung, and F. Buscot. 2004. 
Characterization and spatial distribution of ectomycorrhizas 
colonizing aspen clones released in an experimental field. 
Mycorrhiza 14:295-306. 

Karlinski, L., M. Rudawska, B. Kieliszewska-Rokicka, and T. 
Leski. 2010. Relationship between genotype and soil 
environment during colonization of poplar roots by 
mycorrhizal and endophytic fungi. Mycorrhiza 20: 
315-324. 

Kennedy, P. G., K. G. Peay, and T. D. Bruns. 2009. Root tip 
competition among ectomycorrhizal fungi: are priority 
effects a rule or an exception? Ecology 90:2098-2107. 

Kleemann, F., M. von Fragstein, B. Vornam, A. Miiller, C. 
Leuschner, A. Holzschuh, et al. 2011. Relating ecologically 
important tree traits to associated organisms in full-sib 
aspen families. Eur. J. Forest Res. 130:707-716. 

Labbe, J., V. Jorge, A. Kohler, P. Vion, B. Marcais, C. Bastien, 
et al. 2011. Identification of quantitative trait loci affecting 
ectomycorrhizal symbiosis in an interspecific poplar cross 
and differential expression of genes in ecotmycorrhizas of 
the two parents Populus deltoides and Populus trichocarpa. 
Tree Genet. Genomes 7:617-627. 

Lang, C, J. Seven, and A. Polle. 2011. Host preferences and 
differential contributions of deciduous tree species shape 
mycorrhizal species richness in a mixed Central European 
forest. Mycorrhiza 21:297-308. 

Lapierre, C, B. Pollet, M. Petit-Conil, G. Toval, J. Romero, 
G. Pilate, et al. 1999. Structural alterations of lignins in 
transgenic poplars with depressed cinnamyl alcohol 
dehydrogenase or caffeic acid 

O-methyltransferase activity have an opposite impact on the 
efficiency of industrial kraft pulping. Plant Physiol. 1 19: 
153-163. 

Leple, J. C, A. C. M. Brasileiro, M. F. Michel, F. Delmotte, 
and L. Jouanin. 1992. Transgenic poplars - expression of 
chimeric genes using 4 different constructs. Plant Cell Rep. 
11:137-141. 

Matson, P. A., W. J. Parton, A. G. Power, and M. J. Swift. 

1997. Agricultural intensification and ecosystem properties. 

Science 277:504-509. 
Molina, R., H. Massicotte, and J. M. Trappe. 1992. Specificity 

phenomena in mycorrhizal symbiosis: community-ecological 

consequences and practical applications. Chapman and Hall, 

New York. 

Monk, K. A., and G. J. Samuels. 1990. Mycophagy in 

grasshoppers (Orthoptera, Acrididae) in Indo-Malayan rain 
forests. Biotropica 22:16-21. 



© 2012 The Authors. Published by Blackwell Publishing Ltd. 



1947 



Fungal Communities of a Poplar Plantation 



L. Danielsen et al. 



O'Brien, H. E., J. L. Parrent, J. A. Jackson, J. M. Moncalvo, 
and R. Vilgalys. 2005. Fungal community analysis by large- 
scale sequencing of environmental samples. Appl. Environ. 
Microbiol. 71:5544-5550. 

Pena, R., C. Offermann, J. Simon, P. S. Naumann, A. 
Gessler, J. Hoist, et al. 2010. Girdling affects 
ectomycorrhizal fungal (EMF) diversity and reveals 
functional differences in EMF community composition 
in a beech forest. Appl. Environ. Microbiol. 76: 
1831-1841. 

Perry, B. A., K. Hansen, and D. H. Pfister. 2007. A 
phylogenetic overview of the family Pyronemataceae 
(Ascomycota, Pezizales). Mycol. Res. 111:549-571. 

Peter, M., F. Ayer, S. Egli, and R. Honegger. 2001. Above- and 
below-ground community structure of ectomycorrhizal 
fungi in three Norway spruce (Picea abies) stands in 
Switzerland. Can. J. Bot. 79:1134-1151. 

Pilate, G., E. Guiney, K. Holt, M. Petit-Conil, C. Lapierre, 
J. C. Leple, et al. 2002. Field and pulping performances of 
transgenic trees with altered lignification. Nat. Biotechnol. 
20:607-612. 

Polle, A., and C. Douglas. 2010. The molecular physiology of 
poplars: paving the way for knowledge-based biomass 
production. Plant Biol. 12:239-241. 

R Development Core Team. 2009. R: a language and 
environment for statistical computing, (ed. Computing 
RFfS), R Foundation for Statistical Computing, Vienna. 

Reverchon, F., P. M. del Ortega-Larrocea, and J. Perez- 
Moreno. 2010. Saprophytic fungal communities change in 
diversity and species composition across a volcanic soil 
chronosequence at Sierra del Chichinautzin, Mexico. Ann. 
Microbiol. 60:217-226. 

Rodriguez, R. J., J. F. White, A. E. Arnold, and R. S. Redman. 
2009. Fungal endophytes: diversity and functional roles. 
New Phytol. 182:314-330. 

Schulz, B., and C. Boyle. 2005. The endophytic continuum. 
Mycol. Res. 109:661-686. 

Smith, S. E., and D. J. Read. 2008. Mycorrhizal symbiosis. 
Academic Press, London. 

Stefani, F. O. P., J. M. Moncalvo, A. Seguin, J. A. Berube, 
and R. C. Hamelin. 2009. Impact of an 8-year-old 
transgenic poplar plantation on the ectomycorrhizal 



fungal community. Appl. Environ. Microbiol. 75: 
7527-7536. 

Taberlet, P., L. Gielly, G. Pautou, and J. Bouvet. 1991. 

Universal primers for amplification of 3 noncoding regions 
of chloroplast DNA. Plant Mol. Biol. 17:1105-1109. 

Tagu, D., P. Faivre Rampant, F. Lapeyrie, P. Frey-Klett, P. 
Vion, and M. Villar. 2001. Variation in the ability to form 
ectomycorrhizas in the Fl progeny of an interspecific poplar 
(Populus spp.) cross. Mycorrhiza 10:237-240. 

White, T. J., T. Bruns, S. B. Lee, and J. W. Taylor. 1990. 
Amplification and direct sequencing of fungal 
ribosomal RNA genes for phylogenetics. Academic Press, 
New York. 

Supporting Information 

Additional Supporting Information may be found in the 
online version of this article: 

Figure SI. Scheme of the poplar plantation used for the 
field trial (A) and photo of a typical soil core (B). 
Figure S2. Rarefaction curves of (A) soil, (B) root, (C) 
ECM root tip, and (D) combined samples. 
Figure S3. Detected genera richness per fungal families 
clustered along their presence in soil and root samples. 
Table SI. Relative abundance (mean ± SD) of the 198 
fungal families detected by 454 pyrosequencing analysis in 
soil and root samples. 

Table S2. List of fungal families used for cluster analyses 
(see Fig. 3). 

Table S3. Fungal species were classified into different 
groups according to their ecological lifestyles. Literature 
source is indicated. 

Table S4. Fungal species detected on ECM root tips of 
poplar plants by morphotyping/ITS-sequencing. 

Please note: Wiley- Blackwell are not responsible for the 
content or functionality of any supporting materials sup- 
plied by the authors. Any queries (other than missing 
material) should be directed to the corresponding author 
for the article. 
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